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In large hot tokamaks like JET, the width of the reconnecting layer for resistive modes 
is determined by semi-collisional electron dynamics and is much less than the ion Larmor 
radius. Firstly a dispersion relation valid in this regime is derived which provides a unified 
description of drift-tearing modes, kinetic Alfven waves and the internal kink mode at low 
beta. Tearing mode stability is investigated analytically recovering the stabilising ion orbit 
effect, obtained previously by Cowley et al. (Phys. Fluids 29 3230 1986), which implies large 
values of the tearing mode stability parameter A' are required for instability. Secondly, at 
high beta it is shown that the tearing mode interacts with the kinetic Alfven wave and 
that there is an absolute stabilisation for all A' due to the shielding effects of the electron 
temperature gradients, extending the result of Drake et. al (Phys. Fluids 26 2509 1983) to 
large ion orbits. The nature of the transition between these two limits at finite values of 
beta is then elucidated. The low beta formalism is also relevant to the m — n = 1 tearing 
mode and the dissipative internal kink mode, thus extending the work of Pegoraro et al. 
(Phys. Fluids B 1 364 1989) to a more realistic electron model incorporating temperature 
perturbations, but then the smallness of the dissipative internal kink mode frequency is 
exploited to obtain a new dispersion relation valid at arbitrary beta. A diagram describing 
the stability of both the tearing mode and dissipative internal kink mode, in the space of 
A' and beta, is obtained. The trajectory of the evolution of the current profile during a 
sawtooth period can be plotted in this diagram, providing a model for the triggering of a 
sawtooth crash. 
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I. INTRODUCTION 

Magnetised plasmas are potentially subject to electromagnetic tearing mode instabilities. In 
cylindrical (or slab) geometry these instabilities are driven by the potential energy associated with 
plasma currents flowing along the magnetic field, characterised by a quantity A' and occur 
when magnetic reconnection due to dissipative effects can take place in a narrow region around a 
"resonant surface", i.e. where /cy, the wave number of the tearing mode along the magnetic field 
direction, is zero. Such surfaces can occur in the presence of magnetic shear. The stability of 
this mode is given by a dispersion relation A' = A'^(u}), where A'^{uj) is a function that depends 
on the plasma response at the resonant layer and cj is the mode frequency. This instability is of 
importance in both laboratory plasmas, such as those associated with magnetic confinement fusion, 
and astrophysical plasmas. 

Although many reconnection events involve nonlinear phenomena, two situations in which linear 
theory can be important are sawteeth and mode locking phenomena in tokamaks. In present 
day tokamaks such as JET and the international burning plasma device ITER, currently under 
construction in France, the m = 1, n = 1, tearing (or "reconnecting") mode, and related internal 
kink mode, where m and n are the poloidal and toroidal mode numbers, are believed to trigger the 
sawtooth instability, a periodic disturbance of the core of the toroidal plasma column that occurs 
when the safety factor, q, falls below unity and degrades performance. The rapid sawtooth event 
that is observed in hot tokamaks Q| could be associated with a high threshold in A' [sometimes 
replaced by Xh = — vr/AV^, where Xh is proportional to the growth rate of the ideal MHD internal 
kink mode and q{rs) = 1], followed by some nonlinear destabilisation mechanism that suddenly 
releases the available magnetic energy. For example, the collisional Braginskii two-fluid model 
predicts a stable window in plasma pressure between the drift-tearing and the internal kink modes 
J|, so that the sawtooth crash may be associated with crossing the higher pressure instability 
boundary. The Porcelli et al. model [5| employs a sawtooth trigger criterion given by the conditions 
for the onset of growth of the resistive internal kink mode p, Q|, but it involves some kinetic 
effects (a large Larmor radius model for the ions and a rapid electron transport along the magnetic 
field). The resulting criterion can be expressed in terms of a critical value of the magnetic shear, 



s = rd{lnq)/dr, at the rational 
sawtooth behaviour in ITER 



= 1 surface. Such qualitative models are often used for describing 
8|. Determining the critical value of A' for the instability of these 
modes in such plasmas, as caculated from a relevant theory, should help to model and control the 
sawtooth activity in these devices. 
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For mode locking phenomena, one needs to calculate the low level of magnetic perturbation 
induced by an external coil at a mode resonant surface that produces an electromagnetic torque 
sufficient to overcome viscous plasma forces 0, [l^ . This torque can be calculated in terms of the 
layer quantity A^(a;), but evaluated at a frequency lo representing the mismatch between the E x B 
and diamagnetic rotation of the plasma at the layer and the static external coil, rather than the 
natural frequency of a tearing mode. This information is very important for the resonant magnetic 
perturbations being considered for ELM control on ITER [ll]. 

Initially, tearing modes were studied theoretically using the equations of resistive MHD in cylin- 
drical geometry, when instability occurs if A' > although a later toroidal calculation ob- 



tained A' > A^^^^, where A^^..^^ is a critical value related to the favourable average curvature |12|. 
However, for hotter plasmas various kinetic effects become important. The inclusion of diamag- 
netic effects [13I leads to a reduction of the growth rate and imposes a rotation with frequency 
UJ ~ uj^^e = kyTe/{eBLn), where ky is the wave-number in the magnetic surface but perpendicular 
to the magnetic field (in a tokamak ky = m/r ), Tg is the electron temperature, e is the electron 
charge, B the magnetic field strength and L.„ = [d{\iine) / dr]~^ is the density scale-length. In this 
case the mode is identified as a drift-tearing mode. 

In contemporary fusion devices the "semicoUisional" regime is relevant; this is a situation in which 
the electron collisional transport along the perturbed magnetic field becomes comparable with the 
mode frequency, i.e. uj ~ k'^^vff^^/ue, where vthe = (2T'e/me)^/^ is the electron thermal speed, rrie 
being the electron mass, and is the electron collision frequency. Since A;|| = kyx/Lg, where Lg is 
the magnetic shear length (in a tokamak Lg = Rq/s with R the major radius), and x the distance 
from the resonant surface, this defines a semi-collisional width 60 = Ls{u}^:ei'e)^^'^ / (kyVthe) 
When 5o ^ Pi, where pi = {2miTi)^/'^ /eB is the ion Larmor radius, rrij and Tj being the ion mass 
and temperature, respectively, a small Larmor radius approximation can be employed for the ions, 
but the opposite limit can be more relevant. 

Three key parameters govern such calculations: /3 = Pe{Ls/ Ln)'^ /2, where (3^ = 2pQnf,Te/ B'^ 
is the normalised electron pressure, a coUisionality parameter C = d.hl{ye/^*e){TrLe/^i){Ls/Ln)^ , 
and A'. To see the significance of /3 we note that the collisional resistive layer width, (5^, i.e. the 
width of the current channel at the resonant surface, can be estimated by balancing the mode 
frequency against the resistive diffusion across this width: u ??||/'^^, where 7y|j is the resistivity 
parallel to the magnetic field. For drift-tearing modes with uj ~ uj^e we find that the ratio of 
resistive and semi-collisional widths is given by ((5,j/(5o)^ ~ 2//3, so that /3 is the key parameter 
controlling semi-collisional effects: for /3 ^ 1 the resistive layer width is much narrower than the 
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semi-collisional one so the semi-collisional effects dominate. 

Regarding the parameter C, we can write C = 0.51((^o/Ps)^5 where = T^^'^pi with r = T(,/Ti 
(i.e. pg is the ion Larmor radius based on Tg), so it measures the ratio of the semi-colhsional width 
to the ion Larmor radius. Of course, a cold ion model, i.e. r — )■ oo, decouples the parameter C from 
the ratio of the ion Larmor radius to the semi-collisional layer width and it is then unnecessary to 

n 

account for finite ion orbit effects. Drake et al. [18| considered such a model; in the limit C <^ 1 
they found that the tearing mode is stable until A' reaches a critical value: A^,^^^ ~ P~^^$ for low 
values of /3. For /3 ^ 1, they found that the equilibrium radial electron temperature gradients 
shield the resonant surface from the magnetic perturbation and the tearing mode is stable for all 
A'. The cold ion model with ps > 5o h^-s, since then, been revisited by a number of authors using 
several, more or less ad-hoc, electron models. Thus Grasso et al. {l^ considered a low /3, three field 
model including the effects of ion viscosity, particle diffusion and finite electron compressibility; 
while growth rates were affected by electron diamagnetism, the instability threshold was found 
to be unchanged. This model has been recently re-derived by Fitzpatrick [2^ (in the absence of 
density and temperature gradients), and applied to collisionless and semi-collisional situations. In 



a subsequent paper Grasso et al. 



21[ extended their work to a four field model including ion 



motion along the magnetic field, demonstrating the stabilising influence of ion acoustic waves when 
/3e > A' ps{Ln/ Lg)^^^ , discovered earlier by Bussac et al. 



22|. None of these models, however, can 
reproduce the semi-collisional electron conductivity calculated in Ref. {isl because the authors 
ignored electron thermal effects. However, thermal effects including anomalous electron cross-field 
transport have been investigated numerically in Ref. 
been studied 24 



23] . The effect of toroidal geometry has also 



For a high temperature tokamak, however, the condition pi <C 5q is only valid at very low 
values of magnetic shear s, a situation that might nevertheless arise at the (7 = 1 surface where 
the m = n = 1 tearing mode is resonant. When the width of the reconnection layer becomes 
narrower than the ion Larmor radius a fully gyro- kinetic description of the ions is necessary, as 
Antonsen and Coppi {2^, 



noted 



Hahm and Chen 



27 



28|, Cowley et al. 



29|, Pegoraro and 



Schep [30fl, and Pegoraro, Porcelli and Schep As a result of the non-local character of the finite 
ion orbits, the problem involves the solution of a set of integro-differential equations. Cowley et al. 



29| considered both collisionless and semi-collisional electron models. In the collisionless case they 



recovered the strong stabilising effect from the ion Larmor orbits found by Antonsen and Coppi 



26|, but also found a similar effect in the semi-collisional case. This is due to the non-local ion 



response to electrostatic perturbations. Mathematically, it is essential to keep the small corrections 
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arising from the {kpi)~^ tail at high kpi in F{kpi), the Fourier transform of the ion response with 
respect to distance from the resonant surface. Calculations using a Fade approximation to F{kpi), 
as in part of Refs. Q and 01, would fail to capture the effect of this tail. The calculation by 
Cowley et al. {2^, which was restricted to low /3, found the stabilising effect from large ion orbits 
was characterised by A^j,^+~ pj^ rjl^ \n.{pi/6n), where r/e = d(lnTe)/d(lnne). The semi-collisional 

^- 01,13, 



calculations in Refs. [7||, l27fl, |28| and 



30|, were restricted to the limits of negligible and infinite 



parallel heat conduction, missing important thermal effects. In Ref. the large ion orbit theory 
was applied to the stability of the ideal and dissipative internal kink mode in a tokamak, which 
resulted in higher growth rates than in simpler fluid models |3l| . 



In this paper we revisit the semi-collisional situation in the pi ^ 6q limit addressed in Ref. 
[2J1, but by using a different mathematical approach (based on Fourier transforming the governing 
integro-differential equations) extend it to finite (3 so that we can follow the transition from the 
stabilising effect on tearing modes of finite ion orbits at low (3 to the stabilising effect of the electron 
temperature gradient at high (3. The problem can then be reduced to the solution of a fourth order 
system of differential equations. We take advantage of the separation of scales between the ion 
Larmor radius and the semi-collisional layer width when d^/pi <^ 1, i.e. C <^ 1, to reduce this 
system to the consideration of two separate regions: one the "ion region", where the effects of the 
finite ion Larmor radius appear and another, the "electron region", where the semi-collisional effects 
on the electrons are dominant. The long wavelength limit of the solution in the ion region must 
be matched to an ideal MHD solution characterised by A'. At short wavelengths the ion region 
solution must be matched to the solution in the electron region, leading to a dispersion relation. 

The paper is organised as follows. In Section |Il] we formulate the tearing mode stability problem 
and the separation into the ion and electron regions. In Section IIIII we obtain a unified analytic 
dispersion relation for the coupled drift-tearing mode and the FLR modified internal kink mode in 
the limit /3 <^ 1. When A' > this dispersion relation exhibits coupling between the drift-tearing 
branch and kinetic Alfven waves and we analyse the resulting stability of these modes as a function 
of the parameter A'pif3. In Section [IV] we obtain an analytic solution in the opposite limit, (3^1. 
The transition from instability at low /3 to stability at high /3 is explored in Section |Vl using a finite 
/3 analysis, but one which is strictly valid only for a unique value of r/g. These results on tearing 
mode stability are summarised in Section IVTl Section [VIII analvses the stability of the internal kink 
mode, emphasizing the dissipative mode which is relevant when A' > 0. As well as applying the 
unified low /3 dispersion relation to this mode, we also exploit the smallness of its frequency to 
obtain a dispersion relation valid for finite /3 in the Appendix which is also studied in Section IVIII 
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In Section rVIIII we address the imphcations of the stabihty analysis for sawtooth modelling. Finally 
the results are discussed and summarised in Sec. IIXI 



II. SEMI-COLLISIONAL STABILITY EQUATIONS 

We consider a plasma immersed in a sheared magnetic field in slab geometry, Bq = Bq{z + 
x/Lsy), which serves as an approximation to a cylindrical plasma provided we ignore the effects 
of drifts due to the cylindrical magnetic curvature. Thus Bq is the strength of the guide field and 
Lg the shear length. A gyro-kinetic description of the ions allows us to model finite ion Larmor 
orbits, since we consider these to be large in comparison to the width of the reconnecting layer 
where the electron current flows, pi > 6q; however we do not include the effects of ion motion 
along the magnetic field. Since we do not attempt to describe toroidal geometry, the effects of 
finite ion orbits due to inhomogeneous magnetic fields and their neoclassical consequences cannot 
be addressed. The electrons are described by a semi-collisional fiuid model based on an Ohm's law 
with a thermal force and an electron thermal equation in which parallel electron thermal transport 
competes with the mode frequency. Again, the limitation to slab geometry precludes neoclassical 
effects arising from trapped particles. These ion and electron models allow us to calculate the 
perturbed charge densities and current produced by perturbations in the electrostatic and vector 
potentials and hence, by introducing these in the quasineutrality condition and Ampere's law, we 
arrive at a coupled set of eigenvalue equation. 

Such a model has been provided by Cowley et al. [29| and we reproduce it here: pi > 5 



and 



dpe'P-ip{p)F{ppi). (2) 



dx"^ 5"^ X 

Here the electrostatic potential $ has been scaled to 93 = (uLg/ky)^, A is the parallel component 
of the vector potential, x is the distance from the resonant surface, A:|| = 0, and ip{p) and A{p) are 
the Fourier transform of ip and A. In addition we have uo = 1 — (1 + 1.71rie)/uj, cti = di{l — 1/w), 
uj = ujjijJife-, = kycToe/ieBoa), with a^^ = dlnnoe/dx, do = 5.08, di = 2.13, 

2 ^ ^2 ^ g-if ^^2 B = ^ll 

" ^ 7.9 ,..2 — ^ '^"Oi P o r2 ' 



and 



Fipp^i) = ( i + tVfo - 1) - j^pVf (ro - Ti) ; r„ = e-^/„(pV?/2), (3) 



7 



with pf = 2miTi/{e^Bl). Here, noe and Toe are the equihbrium electron density and temperature, 
respectively. Equations ([T])-(I2D are valid when oo ~ uj^,,, <^ fg, where is the electron collision 
frequency, and at the same time ui ~ ^para//ei^t/ie/^ei- Thus electrons are not assumed isothermal 
and ions are gyro-kinetic. 

For our analysis, it is appropriate to introduce X = x/pi and k = ppi. Then F{k) has the 
following limiting behaviours: 



F{k) [t + {l + r],)/uj]k'^/2 for A; ^ 0, 



(4) 



F(k) F^ + fy for A; ^ oo, 
k 



(5) 



with Foo = — (''" + '^~^); /i = [''" + io~^{l — 7?i/2)]/-^/7r. We solve these equations by asymptotic 
matching between region 1, the ion region, where x ~ pi, i.e. X ~ 1, and region 2, the electron 
region, where x ^ S, i.e. X ~ S/pi, assuming pi ^ 6 so that these two regions are widely separated 
in X. As noted in the Introduction, the high— A; tail accounted for by fi/k plays a key role in this 



matching and affects the mode stability 
Following the work of Refs. js^ and 
where 



J 



7|, it is convenient to calculate in terms of the current J, 



dx^ 



, or J = k'^A, 



so that 



i 5 dJ 
Cj'^fiF{k) Pi dk ' 



In region 1 we obtain 



d_ 

Ik 



G{k) dJ 



+ ^2^^^ = 0; G{k) = F{k)-^-^. 
di fc^ di 



F{k) dk 

Equation ([5]) shows that as A; — t- oo, F — t- F^o = —{t + so that 

J{k) ~ a+A;i/2+M + a^k^/'^-^' for A; ^ oo. 



(6) 



(7) 



(8) 



(9) 



where p is defined by 



1 
4 



l + T \Cj 



+ T 1 



(10) 
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To determine a^/a- we must solve Eq. ([8|) and impose the boundary condition arising from the 
matching to the ideal region at small k (corresponding to large x/pi), given in Ref. js^: 



j(A;)^l + _ZL_a;2(r + i±^)A;3^ A: ^ 0. 



(11) 



3A'pi V ^ 

In region 2 we define t = ex.p{—i7r / A){6oU}^^^ / pi)k. We shall find that it is easier to calculate in 
X space using the scaled variable s = e'x.p{i7T /4){pi/6QU}^^'^)X , and then transfrom to t. Firstly we 
obtain 



(12) 



leading, on integrating once, to 



1 d' 

T + dt'^ 



We back-transform this equation to obtain 

d'^A 



ds"^ 



+ a{s^)Cj^l5A = 0, 



where 



(To + (TlS^ 



(13) 



(14) 



(15) 



(16) 



+ dos2 + dis^' 

(io(w) = do + o"o/(t + ^^^) and di(a)) = di + o"i/(r + lo"^) = di{l + t)/{t + Cj~^). An arbitrary 
constant of integration in Eq. would lead to an unphysical 5— function at s = in Eq. (jlSp . 
From Eq. (jlSp we write the corresponding equation for J(s) 



+ cj2/3 J = 



(17) 



that must be solved with the condition that J{s) is finite and even in s at s = 0. For large argument, 
we find 



J{s) ~ 6+5-3/2+^^ + b.s-^/^-^" for s ^ oo. 



(18) 



so that for small argument in Fourier space 



J{t) ~ + c_ti/2-M for t 0, 



(19) 
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where it can be shown that, recahing J(s) is even. 

c_ _ 6+r(-i + ^) 



tan 



TT /I 



which can be matched to the region 1 solution to yield 



2 V2 



a4 



-i7r/4 



2^ 



Pi 



(20) 



(21) 



Thus, to determine c_ and c+, we see that we must solve Eq. (I17p to obtain 6„ and 6+, deduce 
c_/c+ from Eq. (j20p which can then be used in the matching condition (|2ip to provide a dispersion 
relation. This analytic process has fully accounted for the small parameter Jo/Pt- In general, one 
needs to solve Eq. ([8]) with boundary condition (llip to determine a_ and a+, then solve Eq. (flTl) 
with the appropriate boundary condition at s = to determine 6+ and and hence obtain c_ 
and c+ using result (j20p . Normally this would require numerical solutions of Eqs. ([8]) and (jl7p . 
involving the parameters /3, r, rje, rji and pjA'. However here we shall make some analytic progress 
in the two limits /3 ^ 1 and /3 ^ 1, while also obtaining results valid for finite /3 based on the fact 
that cj — )• 1 for the tearing mode and ^ 1 for the dissipative internal kink mode at large A'/9j. 



III. THE LIMIT p < 1 



A. A unified dispersion relation at low /3. 

Following Ref. js^, a solution of Eq. ([8]) for region 1 satisfying condition ([TT]) at low k can be 
found by expanding in 13 to obtain 



j{k) = exp 



■^2^^ du F{u) 
di Jq u G{u) 



di A' Pi Jo G{u) 
di ) Jo G[u) Jo v^G[v) 



At large k this takes the form 



+ 



vr /I 



7 



^^'(I + t) 

1 + OJT 



A'pi \A ^ J 1 + TCdT + l 
f \ X ^ ( 1 I 7"') 

/(ci, T, r/i) - ( - - ptM — — ^ — /(w, T, rii)k for 



where 



I{uj,T,r]i 



dk 



F{k) 1+uT {Co - + COT - r]i/2) 



+ 



G{k) 6j{1 + t) ^uj^{l + T)^{l + k) 



oo, 



(22) 



(23) 



(24) 
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where we have recalled the asymptotic form of F{k) at large k given in Eq. ([5]). We note that the 
term involving fi gives rise to the In A; term in Eq. (I23p . 

In region 2 we solve Eq. (I17p to zeroth order in /3 to obtain 



(25) 



where Jq is a constant. In next order we find 



(26) 



At large s, evaluating the integral = a (s^) by contour integration, we obtain 



J{s) ~ Jo 



r (7l S + W in S 

2di(s+ + s-j V'5+'S- / di 



for s — )■ oo. 



(27) 



where 



Oq 



2di ^ Y ^1 



At low /3 we must take the limit /-t — )■ 1/2 of J{s) in Eq. (jl8p : 



1 



J(s)~^|6_ (^l + w^/3^1ns^ + ...j +6+(s + ...)[>, fors^oo, 

so that we can identify and 6_ by comparison with Eq. (j27p : 

^ vr 1 /do 



2(ii 5+ + S- 
Then Eq. (|2ip in the limit ^ 1/2 implies 

c_ 2di 1 b+ 



0-1 



so that 



(28) 



(29) 



(30) 



(31) 



Jit) w c_t^/2"^ + c+t, for t ^ 0. 



(32) 



Matching the powers in the solutions given by Eqs. (|23p and (j32p . with the use of Eq. ()3ip . 
yields the dispersion relation 



e 4 . 



(5o w^/^cJi \ldQ + 2\/di 



Pi y di a^y di + ai 



1+Cjt , 

1 + T ' 



1) ^-ci(a;-l)/3^/ 



(33) 



1 - w(cu - 1) 



7r/3 



(cj-l)(i:Dr+l->?»/2) 



ln(e*4^ 
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To analyse Eq. (|33p it is useful to cast it in the form 



5 —A{uj)B{6j) + 6jVI + 6jtC{uj)D{uj) = 0, (34) 



where we have defined 



A^uj) = Vl + T\J w - (1 + 1.7l7?e) + doil + Cot) + 2v^di(l + r)(l + a)T)(j 

C{uj) = 1 ^^^uj{Cj — T, r/j) 

vr 

= [uj-{l + 1.71r/e)]\/w(l + T) + (w - l)v/(ii(l + i2'T). 
The important ln(/9i/5o) term in B[Cd) arises from the fi/k correction in Eq. Since 5q/ pi <^ 1, 
there appear to be four uncoupled solution branches for the lowest order eigenvalue, ujq, given by 

LjQ = -1/r, C(wo) = 0, and D{ljo) = 0. (36) 

The first of these is an ion drift mode. The fourth branch, D{ujq) = 0, is the drift-tearing mode 
whose stability is discussed shortly. The other two branches arise from solutions of C(a)o) = 0, 
corresponding to a pair of kinetic Alfven waves (KAW's) if A' > 0. To see this we consider the 
limit \uj\ ^ 1, when ~ 0.798 + (1.194 + 0.399rji)/u}, where we have taken the limit r = 1 for 
simplicity. Then we find that 



cd = ±J (37) 

where we require /3~^ ^ ^'Pi ^ 1, so that the conditions o)^ <C 1 and cD ^ 1 are met. In Alfvenic 
units, the frequency of the mode is 



a; = ±a;^^M. = ±^^;^V.g, (38) 

where we have defined the Alfven speed va = B / ^/rnjul. The two solutions in Eq. (j38p correspond 
to one KAW propagating in the electron diamagnetic direction and another in the ion direction. 

Including the 6q/ pi correction to Eq. (134^ . we find that these two KAWs are stable. The lowest 
order drift-tearing and KAW mode frequencies as functions of A'piP are shown schematically in Fig. 
[TJ At large values of A', a further solution of C(w) = is associated with the dissipative internal 
kink mode discussed in Section [VIH while negative values of A' give rise to the FLR modified ideal 
internal kink mode analysed in Ref. jsol . 

However we note that at small /3 when 5o/pi ~ /3^, but A'/9j/3 ~ C(l), which introduces an 
interaction between the drift-tearing mode and the KAWs, all four branches become coupled and 
Eq. ([34|) no longer has a natural expansion parameter. 
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Figure 1: The solutions of the dispersion equation for the drift-tearing mode and two KAWs as a 
function of A'pil3/n. There is also an ion drift mode. The drift-tearing mode and the KAW propagating in 
the electron direction are strongly coupled in the vicinity of A'pif3 = {A' pi(3)c given in Eq. ([5^ . 

B. The drift-tearing mode 



Returning to the drift-tearing mode and assuming A'pif3 <C 1 so that we have an uncoupled 
mode, the lowest order solution, D{ijj) = 0, is given by 




Figure 2: The drift-tearing mode frequency [the solution of Eq. (|39p ] as a function of ije for r = 
negative r]e < —1.17 a second solution exists. Also shown as a dotted line is the fit 



1. For 
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given by 



iu = luq ^ 1 + 0.7Ar]e. 



(40) 



Figure [2] also shows a second solution for negative r]e, i.e. reversed density gradients or hollow 
temperature profiles; in fact it exists for 7]^ < —(1 + t^^)/1.71, i.e. r]e < —1.17 for r = 1, and for 
larger negative rj^ it can also be fitted by the form (|4U|) . Thus, in the limit of a flat density profile, 
one does obtain the same frequency, irrespective of the sign of rjg. The growth rate follows from 
next order in do/pi : 



X -1/2 

„ _ 7 OquJq , , , . 

^=7, ~ ~ ^^r-it) 

tO*e ~ 



7r/3 



where A[,^^j is the critical A' for instability 



^crit = -2n J. \2 r'o^ + 1 - ^ K - 1) In e 

Pi uj^{l + r)^ V 2 / \ 



(41) 



^{000 - 1)1, (42) 
Pi 



provided r/j/2 < ujqt + 1, where we note that expression (j40p implies the factor (cDq — 1)^ oc r^g. This 



result closely resembles that obtained in Ref. 



29|. The integral I[ujQ{r]e),r]i,T], which is negative. 



scales approximately as I 
stabilisation found in Ref. 



1/2 

rje , so the stabilising effect of this term resembles the diamagnetic 



A'i^ = ^(cio-l)|/| 



Pi 



(43) 



as can be seen by taking the r — )• 00 limit of Eq. (I43p . noting that I ^ t ^^'^ in this situation. The 
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Figure 3: The scaling of the quantity A^^'^" in Eq. ([43] 



quantity A^'J^jj is shown as a function of r/g in Fig. [3l The two different effects are shown in Fig. HI 
where 
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Figure 4: The diamagnetic and ion- kinetic contributions to the critical A' [given in Eq. (HH)) and (|44p . 
respectively] as a function of the electron temperature gradients 7]^. Here /3 = 10^^, (5o/Pi = 10^'^, r = 1, 

It should also be noted that, in principle, a large enough ion temperature gradient {7]i/2 > (2}ot+ 1) 
can give a negative critical delta prime, however this requires values of rji that are probably only 
physically relevant in the flat density limit. 

When < we have seen there are solutions of Eq. (|39p with cD < 1, in which case the integrand 
over k in the integral /(a)) has a pole at A; = ko, where G(fco) = 0. This zero in G{k) corresponds 
to a resonance with an electrostatic drift wave of wave-number k. The integral I{u}) then has an 
imaginary contribution: by considering co to have a small imaginary part, we find the integral 
along real k must pass above the pole at = /cq to ensure causality. This contribution provides 
a stabilising effect that increases the value of A^^-^ in Eq. (j42|l via damping on the electrostatic 
electron drift wave. 



The kinetic Alfven wave. 



According to Eq. ()4ip . as the tearing parameter A'pi/{TT/3) is increased beyond its marginal 
stability value (j42p . the drift-tearing mode becomes increasingly unstable. However, when A'pif3 ~ 
0{1) the approximation C{uJo) ~ 1 fails; the dispersion relation can then be written in the form 



e 4 — —A[oo) hwyl + wr 



TT 



1 



TT 



Ld{uj — l)/(w) 



D{ib) = 0. 



(45) 
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At the drift-tearing frequency, cDq ~ 1 + 0.74r/e, the ion integral I is real and positive so the KAW 
propagating in the electron direction becomes strongly coupled to the drift-tearing branch when 
A'piP ~ At yet higher values of this parameter the modes again decouple. We shall show 

below that the KAW propagating in the electron drift direction connects with the tearing mode 
having the frequency ojq. Because C{ujo) has now changed sign, the branch with u = loq has become 
stable. On the other hand the drift-tearing mode couples to the KAW given by the solution of 
C{ujq) = 0. Since one can show that, near a) = 1, 



(46) 



this branch has an unstable solution with 



1 + 



8V2y^do + 2Vd^ - 1.71r?e/(l + r) 1 6o 



1.7lrje{l + T + r]i) 



/3^A'pi Pi 



(47) 



and a growth rate 



7 



1 



16 <5o Vdo + 2VdI - 1.7l7?e/(l + r) 
V27r Pi /33AVa.71r?e(l + r + 77,) 

1 6o A' Pi Vdo + 2VdI - 1.71r7e/(l + r) 



(48) 



Thus, as P increases, w — )■ 1 and the growth rate decreases to zero. In fact, at very high values of 
A' Pi the mode becomes stable. The condition for this is 



A' Pi > 2.427r 



Pi 



(49) 



'^0 Vdo + 2V^- 1.7l7?e/(l + r)' 

This stabilisation mechanism does not exist when /3 < (So/pi)^^^ since the perturbative solution 
of Eq. (|45p fails for A' pif3 ~ 0{1). Numerical solution shows that the unstable drift-tearing mode 
persists as A' pi increases beyond the critical value given in criterion (|42p . eventually merging with 
the dissipative kink instability discussed later in Section [VIII when A' pi S> 1. 
We also note that Eq. (j^H]) implies stability when 



do + 2Vdr- 1.71r7e/(l + r) = 0, 



(50) 



corresponding to a critical rje- Hence, as w — )• 1, the drift tearing mode is also stable if 



r?e > 4.68(1 -hr). 



(51) 
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D. The coupled modes 



To analyse how the two branches interact we examine the solution of the dispersion relation (I45p 
in the vicinity of the cross-over point. At this point, C{uj) = D{Cj) = 0, so that a) is given by the 
solution of Eq. ()39p . i.e. Cj = cjo, and the critical value of the parameter IS.' piP is given by 

TT 



(A'p./3), 



(52) 



d)o(a)o - l)^('^o,'?/i,T)' 

where the integral / is real and positive for modes propagating in the electron direction. Here the 
two modes interact strongly. After introducting 6uj = Co — ujq and 5{A' pif3) = A'pi(3 — {A'pif3)c, a 
perturbation theory yields 



dC dP 

dib duj 



5(AV^ _,.8oA{Cjo)B{Cjo) 

; he 4 . 



Pi woa/I + '^gt 



0. 



(53) 



One can trace how the branches reconnect about this point using the solution of Eq. (i53l) : 



5Cj± 



1 



5{A'pS) 



2{dC/du;h, {A'pS)c 2(9C7/5(i)i, 



When 



S{A'p,/3) 
{A'p,p)c 



Pi {dD/dLb)c,„ u:Q^Jl + LOQT 

2r(5o (dC/du;)^,^ A{ujo)B{u}o 



(54) 



4e" 



solutions (I54p simplify to 



6uj± 



5{A'piP) 



Pi {dD/du)cjo cjoVl + t^oT 



5{A'p,^) 



(55) 



=Fe" 



2{dC/du)^, {A'pS)c 
5q A{i1jq)B{i1jo) 



2{dC/du^)c,, {A'pS)c 



Pi woVl + 



5{A'pS) 



(56) 



sign [{dC/dQj)c, 



2idD/du;)^, {A'pS)c 

Thus we see that for A' pij3 < {A'pif3)c, 6(2)+ connects to the stable KAW branch, while for A'pif3 > 
{A' pi(3)c it connects to the stable drift-tearing mode, whereas 5a)_ connects to the unstable drift- 
tearing mode for A' pij3 < (A'pj/3)c and to the unstable KAW for A' pi(3 > (A' pi/3)c. This is shown 
in Fig. [5l where we plot 5(2)± in the vicinity of the cross-over, by varying the parameter A'piP 
at constant /3 = 0.2, for 60/ pi = 0.1, and r]e = r]i = t = 1. [the quantities in Eq. ()56p are very 
insensitive to these latter three parameters]. Very close to A' pi(3 = {A'pi(3)c the unstable branch 
has a larger growth rate than is given by Eq. (j48p . namely 
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(y-y)/y= 



Figure 5: The variation of the real frequencies of the two modes with pi f3, for fixed /? ~ 0.2 and / pi = 0.1, 
in the vicinity of the cross-over point, y = A'pil3 — {A'piP)c = Uc and uj ~ ujo, showing how the KAW 
converts to the drift-tearing mode and vice versa. These results are for a typical case: t],, ^ rji = t = 1. 
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0.2 



(y-y.)/yc 



0.0 




Figure 6: The equivalent variation of the growth/damping rates of the two modes in Fig. [Sj showing how 
they peak at the cross-over point. These results are for rje ~ iji = t = I. 



In Fig. [6] we show the dependence of the growth (or damping) rates of the two modes for the same 
parameters as in Fig. [5l 



IV. THE LIMIT /§ > 1 

The analysis in the previous section was limited to the case /3 <C 1; in this Section we explore 
the opposite case: /3 ^ 1. We anticipate that the "intermediate eigenvalue" fi will be finite, so that 
at high /3 we will have a) — 1 ~ 0(/3~"^), as in Ref. [isl . However, we have also seen in the low /3 
theory of Section ITTIl that a) — )■ 1 when A'pif3 ^ 1. In Section |V] we will consider the fact that we 
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can also have a)— l~e^lat finite /? to provide a bridge between the high and low /3 limits. 

If w — 1 ~ 0{j3~^), then a/diF{k) ^ 1 in Eq. ([8]) for region 1, unless A; ^ 1, in which case 
F{k) oc /c^.Thus for finite k [not just for A; 3> 1, as in Eq. ([9])] the solution of Eq. ([8]) is 



(58) 



For A; ^ 1, however, Eq. ([8]) takes the form 



dk 



2{1 - u;-^) + {1 + r]i + T)k'^ 
(1 +r/i + t)A;2 



dJ 



dk 



^+ 7-/^' To =0. 



J 



A;2 



(59) 



This has a solution 



J(fc) = ai2Fi(-i + |,-i-^;-i;y 



where 2-^1 is the hypergeometric function |32| . 



(60) 



y 



1 + ??i + r 



(61) 



2(1 -w-i)' 

and «! and 02 are two constants to be determined by matching solution (160p to the ideal MHD 
boundary condition (jlip for A: ^ 1. This matching condition yields 



02 _ ^ 

ai 3k A'pj 
The large k behaviour of solution (|60p is given by 



(62) 



01 



r(-i/2)r(-^) 
r2(-i/4-M/2) 



ia2 



+ y 



l/4+Ai/2 



T2(-l/4 + /i/2) 

where T is the Euler Gamma function 
the quantites d± in Eq. ([9]) are related by 



r(-i/2)r(/.) . r(5/2)r(M) 

oi— 757 — —7- 7— — ^a2 



r(5/2)r(-^) - 

r2(5/4-^V2). 

for A; — )• 00, 



T2(5/4 + ///2) 



(63) 



33|. Thus, the asymptotic form given by Eq. (|63p implies 



a_ 

0.4 



-2fi 



r^(-M/2-l/4) ^8ai r^(-At/2+5/4) 



(64) 



r^li/2-1/4) '8ai r^(/./2+5/4) 

In region 2 we follow the treatment of Ref. (isl to solve Eq. (I17p . obtaining a solution that is small 
at s = : 

C.2 



J{s) 



(To + 



l + dos^ + dis^V Si ^Vs 



(65) 
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where K^^ is a modified Bessel function 



J{s) 



and = l.TlT^e/S/di. For large s this has the form 

3/2+M -j^ 



T{l-fi)\stJ T{l + fx)\stJ 

Calculating the Fourier transform of expression (j66p we obtain 

6_ /sA2M r(l-/i) 

2) r(i + /i)' 



25 \ -3/2-/. 



00. 



Thus, with the aid of Eq. (j20p . we have 



c_ 



2V^ r(/.-i/2)r(i + ^) ^ 

— ; — ^ r tan 



^sty r(-/i- i/2)r(i -^) 

Finally, matching region 1 and region 2 solutions, we obtain the dispersion relation 

{^, + l/2)T\-^,) fP-cot [vr(i + f)]l 



where 



and 



(-/i+l/2)r2(;.) I^P-cot [vr(i-f)] 



R 



V2 (1 + r 

vr (a; - Iji/^ 



1.71r?e (^g 1/4 - ' 
f 7?,)'/' r(5/4 - ^/2)r(5/4 + ^/2) 



r(3/4 - /u/2)r(3/4 + /i/2) ' 
Introducing the collisionality parameter C of Ref. |l8l |. 



C 



0.51-^^^ = 2.04^, 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



equation (jGOp can be recognised as essentially identical to their semi-collisional, high /3, cold ion 
result, generalised to finite Tj/Tg and rji. While surprising at first sight, we see that the solution 
(j60p of Eq. ([8]) only utilised the low k expansion of F(k), corresponding to a small ion Larmor 
radius approximation, although the solution remained valid for finite k. Furthermore this is the 
solution that we would have obtained using a Fade approximation to F{k), as in Ref. [7|; however, 
here it is exact for /3 ^ 1. As is evident from our low f3 result, this approximation fails to capture 
the "1/A; tail" of F{k) at high k in Eq. 1^ that is essential to the result of Ref. 

As shown in Ref. jisl . the dispersion relation ()69p predicts stability for all values of A'. An 
analysis at large A' shows the damping rate becomes small 



7 



-l/(i?i/2A'p,) 



(73) 



The result ()48p . representing the KAW branch with ~ 1 at lower /3, also has a growth rate that 
tends to zero as {A' pii3^)~^ . In the next Section we present a model that demonstrates how the 
transition to stability occurs at finite /3. 
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V. TRANSITION TO STABILITY AT FINITE /3 



As we have seen, the KAW branch has a frequency approaching w ~ 1 as /3 increases, as also 
does the high (3 theory, and we exploit this property to explore the effect of finite /3, i.e. /3 ~ 0{1). 
For finite s, Eq. (jl7p takes the form 



-^{l + [do- 1.71r/e/(l + t)]s^ + dis^} J - 1.71r/e/3J = 0. 



(74) 



If we consider the special case 



do - 1.71r/e/(l + r) = 2^/dl, (75) 

corresponding to rje = 2.53 for r = 1, and introduce u = d\^^s, then we can write Eq. ([71]) as 

1.71?7e/3 



dv? 



'dis,. 



The solution, even at s = 0, is 



J{s) = (l + u^) '^^^ cos ^\/ 1 — arctan uj , 
which has the large s asymptotic behaviour 



J{s) 



~ s 



d 



tanh (^-Va2 - Ij s ^ s»l. 



(76) 



(77) 



(78) 



When s is such that [Gj — l)s^ ~ 1, the solution for J(s) is obtained in terms of modified Bessel 
functions, as in Eq. (|65p . except that now we no longer impose the vanishing boundary condition 
for small s, but match to the form (1781) instead. Thus we obtain 



J(s) oc 



\^r^^x-{Qj-\)s'^ 

.s7/2 



(79) 



where aiand 02 are constants and we note i^? = 1/4— (a) — l)/3 ~ 1/4. Matching the small argument 
limit of this solution to the form (|78l) . with ~ 1/2 and tl; ~ 1, yields the relationship between 
aiand ■ 



02 _ vr 



1 



A 



tanh (^^\/A2 - 1 



(80) 



At high /3, i.e. high A, we see that 02/^1 0, and we recover the results in the previous Section. 

When (w — 1)5^ ^ 1, the solution (|79l) has the correct form to match to the solution (I63p in 
region 1. As a result we obtain a dispersion relation 

(/X + l/2)r2(-;,) JP-COt [f (i+/^)]l 



2 /^i?/^ = A (/i. A) 



(-/i + l/2)r2(^) 1 P-cot [f (i-/i)] 



(81) 
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which only differs from expression (|69p by the key "shielding factor" A (/i, A) : 



A (yU, A) = 1 — sin(7r/i) 



l-^^^tanhf^y/A^ 
A \ 2 



(82) 



For consistency we need to consider A (/i, A) in the limit n — t- 1/2. We note that then A ~ 1 for 
A > 1, A = for A = 1 and A -4/{ttX^) for A < 1. 

Recalling R ^ 1, we observe that Eq. ([8T]) has a solution with fi ~ 1/2, as required if our 
ordering is to be consistent. This solution is obtained by balancing the large left-hand-side of Eq. 
(jSip by approaching the zero of the denominator on its right-hand-side. Writing /x « 1/2 — /36uj, 
we find 

2 f 2 y ^ 2V2A (1.71 6q 



w ~ 1 + + —rie —, (83) 

l + T + ViKpA'piJ /33/2AV,(l + r + 7?i) V '^1 J Pi 

,»_4^A^(i^!!.)"^_^ (S4) 

Pi\di /3j /3AV,(l + r + 7?,) V 2V2 P^\ di ' j 

which shows that we require I31S.' pi ^ 1 for consistency. In the small (3 limit the growth rate reduces 
to 

^^J2^5q^ 1 / ^/2d}/'AV, Jo\ 

7r^/2Pi 1.71r/e/33AV,(l + r + r?i) V ^-^^ r|Jp^)' 

Taking account of the condition ([75]) on r/g, we see that this coincides with the low (3 result (l48|l of 
the KAW; similarly, the frequency Co reduces to the result (I47p . 

Equation (j84p implies stability if A > 0. This stability boundary corresponds to A = 1, so that 
there is stability for values of j3 satisfying 

/3 > /3, = = 0.34, (86) 

1.71r?e 

for r = 1. Thus the model correctly predicts stability at high j3 as found in Section ITVl If A < 
then there is still stability if A' is sufficiently large: 



A'5oA/^/3r?e>-^. (87) 



For /3 — )• /3c, one finds stability if 



A'5o > ^a/^^. (88) 

4^2 V /3c-/3 
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As P reduces, —A increases and the necessary value of A' for stability decreases. At low f3 the result 
is consistent with the stability condition (l49l) . We note that a similar analysis can be performed 
for another special case: 

do - 1.71r7e/(l + r) = -2y^. (89) 

Now, however, we find the analogue of the screening factor A is always positive so the mode is 
always stable. Condition (j89p . which guarantees stability at finite /3, coincides with condition (j75p . 
the stability condition at low /?, so we have a consistent picture. 

Although our model equation is strictly valid for a specific, though reasonable, value of r/g, it 
captures in a clear way the onset of screening by plasma gradients that stabilises modes in high /3 
plasma. 

VI. TEARING MODE STABILITY DIAGRAM 

In this section we draw together results from previous sections concerning the stability of the 
drift-tearing and kinetic Alfven waves, summarising them in Fig. [71 a diagram in the space of /3 
and A' Pi which delineates stable and unstable regions and regions of validity of these results. 



. 8 
. 6 
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0.2 



10 20 40 50 60 

A Pi 

Figure 7: Stability diagram for the drift-tearing/kiiietic Alfven mode in terms of A'pi and /3. Regions I, 
IV and V are stable. Regions II and III unstable. The finite /3 analysis is not valid in the region below 
the hyperbola $A'pi = 1 (dashed line); the low /3 theory fails below the dotted line /? — {So/piY^^ when 
/3A'P^ - 0{1). Here So/p. = IQ-^. 

The low /3 theory of Section [III B\ valid for f3 <^ 1, shows that the drift-tearing mode is stable 
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in Region I, becoming unstable in Region II when A'/Oj exceeds the critical value given in criterion 
(I42p . When the approximate boundary A'/9j/3 ~ 0{1) is exceeded as in Region III, the analysis in 
Section flll Cl is valid and the kinetic Alfven wave becomes the unstable branch; this remains unstable 
until the criterion (|49p is satisfied and one enters the stable Region IV. Above a critical value of 
/3 = /3c ~ C(l)) given precisely by condition (I86p in Section |V] for the special case of r]e specified 
in Eq. (l75p . the kinetic Alfven wave is stable, as represented by Region V. This stability region 
extends to high f3 as shown in Section IV. The boundary between Regions IV and III asymptotes 
to the critical value of /3 following the scaling A' oc (/3c — /3)~^. The low /3 results discussed above 
are limited by the condition (3 > [So/pi)^^"^ (dotted line in Fig. [7|) when IS.' ~ C'(l). Below 
the dotted line, there is still an instability when A'/jj exceeds the critical value from the criterion 
(j42p . beyond which the growth rate continues to increase with A'pj, eventually merging with the 
dissipative kink mode, discussed in the next Section, when A'pj 3> 1. 



VII. THE INTERNAL KINK MODE 



A. Low /3 theory 



In this Section we analyse the internal kink mode. When A' < we first follow Pegoraro et 
al. Q in introducing the growth rate of the ideal MHD internal kink mode: Xh = —n/A'rs, 
where q{rs) = 1- The solution of the branch C{uj) = of the general dispersion relation (p4|) 
then reproduces the FLR modified internal kink mode dispersion relation obtained, and analysed 
numerically, in Ref. The small terms on the left-hand-side of Eq. (j34p can be retained to obtain 
a correction to the ideal result. However, in the limit Xh — ?• these terms must be retained to 
obtain the dispersion relation for the dissipative internal kink mode which is associated with the 
tD ~ branch of Eq. (|36|) . 

For small values of |tD| and ^ 1, we can approximate the terms in Eq. ()34p : 



Aicj) = VlT^Vdo - (1 + 1.71r/e) = V(1+^K^ 

(90) 

C(w) = l + — ^/3w/(w,r,r7i) 



D{Cj) = -Vdi, 

where now we require pi\Cj\ / {t: (3) ^ 1. Assuming Cj to be real and negative, i.e. in the ion 
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direction, we can calculate the integrals analytically: 



1 - r?^/2 



In 



1 - r?^/2 
-w(l + r)0F 



where 



1 f^dk 



(91) 



(92) 



Go l-r?i/2(l + fe) 

The quantity (1 — rji/2)l2{i]i,T) is plotted in Fig. [8] as a function of rji for r = 1. For consistency. 



(l-r7i/2)l2(r?i,r = l) 



1.25 1.5 



Figure 8: The function (1 — r]i/2)l2 appearing in Eq. ([92|) as a function of rji for t = 1. 



the present analysis of Eq. (j34p is limited to the situation where 6o/pi(3 ^ 1 to ensure |a;| <^ 1. 

Substituting the result (|90p into Eq. (j45p and expressing A' in terms of Xh, we obtain a 
dispersion relation 



a) < In 



(1 + t)V¥ 
1 - ?/i/2 



a; 



2 ; ' V 2 ; Pi 



l±^.(.)(l-?l^e-^l. 



(93) 



2^ 



which is valid for |A///pj| ^ 1. For the special case Xh = this produces, as in Ref. js], an unstable 
mode rotating in the ion direction 

1 So _,-2r 



UJ 



-e 4 . 



(94) 



27rdi V- 2 y ln(p,/3V<5o) p.^^ 
where we have iterated on the ln{—Lb~^) term, ignoring 0(1) corrections. When So/pi <C 
P \Xh \ / Pi we obtain a stable mode in the electron direction: 



UJ 



vrV 2) p,^ln{p,^/\XH\) 



ITT 



ln(p,/3/|AH|) 



(95) 
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Thus there is a marginally stable value of Xh, which can be found by equating real and imaginary 
parts, with cD real, in Eq. (I93p . to yield 



\ crit 



1 



^2 p. 



TT 



2di 



So (1 + r)3/2 y 

The condition that \uj\ remains small at these value of Xjj requires 



+ 7r+(l-^)h{v^)}. (96) 



(97) 



We note that the numerical solution at high A'pi in the region 5o/p« ^ H described in Section 
IIII CI can be followed through the boundary /3 ~ ((^o/Pi)^^^ to merge with the result We 
identify also a strong stabilising effect from rje ■ Eq. (j90p implies the factor d{rie) becomes negative, 
corresponding to stability, if 



r?e > {do - 1)/1.71 ^ 2.38. 



(98) 



Using Eq. (1721) we can express the results ([M|) , (|95p and (|96p in terms of the collisionality parameter 
C of Ref. hjl. 



B. Finite /3 theory 

Because the dissipative internal kink mode corresponds to \uj\ <^ 1 we can seek a solution of Eqs. 
dS]) and ()17p based on an expansion in lj at finite /3, rather that in /3 at finite a), as in Section UTTl 
The details are presented in the Appendix but we note the essential points here. For the ions we can 
exploit the same expansion as in Eq. (|22p . but now based on |a;| <^ 1; in this region of k, F{k) and 
G{k) take simpler forms. However, this expansion fails at large k ~ cb^^ when the finite a) correction 
to G{k) becomes important. The ion response must be calculated in this intermediate region with 
k ~ ib^^ before the solution can be matched to the electron region. The electron solution, i.e. the 
solution of Eq. (I17p calculated in real space, also exploits the condition ^ 1 by considering two 
sub-regions. In the first of these we consider s ~ 0{1) when we take advantage of the fact that the 
quadratic term in the demonimator of the conductivity in Eq. (llSh can be ignored. In the second 
sub-region we take s ~ |(i|~^/2 ^ 1 which again simplifies the equation so that an analytic solution 
can be obtained. We match the solutions in the two sub-regions at intermediate values of s, deduce 
the power-like asymptotic form for s ^ |a)|~^/2 and then use Eq. (|20l) to obtain the solution in 
A;— space. Finally we match this electron region solution to the intermediate ion region solution to 
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provide a dispersion relation for the dissipative internal kink mode. This relation is given by Eq. 
(jA26p in the Appendix, which we reproduce here: 

8^iop,^ di (1 + t)3/2 
We recall that we can express do/ Pi in terms of the collisionality parameter C, [note do/i^Pi) = 



^ 1/2 -/ii r2(Mi) nu;)+ij^e 
1/2 + /.I r2(-^0s(a;) + i4^e- 



(99) 



cosfvr/ii] 



vr 



1 l-r]i/2XH 1 Al-^?^/2 



\/7r 1 + T Pi CO 



1 + r 



-/2 



+ - 



1 + r 



3 , ^, / ^(1+r) . 

Kn + m ^ W 

2 V l-?^/2 



(100) 



with 2fii 
function 



1 + 4/3/(1 + t), ko = V'(l)+^/;(3)-V'(3/2-/ii)-^/;(3/2 + ^i), where ^ is the digamma 



33| . We can solve this analytically in the low (3 limit. Assuming the logarithmic term 
in Eq. (jlOOp dominates, we recover the result (j93p . However Eq. ([99|) allows us to explore the 
effect of finite (3 on the stability of the dissipative internal kink mode. For low values of /3/(l + r) 



Unstable 



Stable 



/3<(5o/Pi)- 

-0.15 -0.1 

Ah/Pi 



Figure 9: The boundary of marginal stability for the dissipative internal kink mode when [Xn/pil ^ 1. The 
complete line is from the numerical solution of the dispersion relation (|99|) . The dashed line is the analytical 
boundary of Eq. (jlOip . the dotted line is the low /3 analytical limit of Eq. ((96|) . Here 60/ pi = 10^^, 
r]e = T = 1, rji = 0. In the region /3 < {do/ Pi)^^^ , the finite /3 theory breaks down. 

there is a threshold value for |Ai^| for the onset of the instability, as in Eq. (I96p . but this is lost 
as /3/(l + r) increases beyond 0.194. It is again possible to calculate the marginal value of \Xh\ by 
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Figure 10: Boundary of the marginal stability region for the dissipative internal kink mode from the numer- 
ical solution of Eq. Solid line: 5o/pi = 10^'^, dashed line: 6o/pi = 10^^. 



equating real and imaginary parts of Eq. (|99|) with real uj (for a consistent solution we must take 
(jj = exp(— 27rz)cDi, with coi real and positive): 



yi + T / d(r?e) 
8 V di 



/3 



vr / sin(7r^i/2) 



1 + T cos(7r;Ui) y sin(37r/ii/2) 



/3 



1 + r 



while a) is given by 



a) = ~H{ni) 



Here 



1/2 + ^1 r2(-^ii) 



/ di 1 - r?j./2 

TTdiVe) (1 + r)3/2 ■ 



cos(7r/.ii) 



1 

2ni 



1/2 -/ii r2(^i) cos(7r/ii/2) ± Y^sin(7r^i/2)sin(37r/ii/2) 



(101) 



(102) 



(103) 



Expression (jlOip correctly reduces to the result (f96|) in the low /3 limit. 

Figure [9] shows the resulting critical value of \Xh \ / Pi as a function of 13 for 5o/pi = 10~^ (other 
parameters are ?7e = li f?i = and r = 1); we recall the low /? limitation to ensure \Cj\ <^ 1 which is 
marked in the figure. The validity regime for the diagram is {Xn/pil <^ 1, i.e. in the region where 
the mode is not ideal, but close to the ideal stability boundary Xn/pi = 0. There is no marginal 
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value of \Xh\ when /3/(l + r) > 7/36 ~ 0.194, corresponding to fii = 2/3 when Eq. (|10ip imphes 
\Xh\ — 5" oo. Figure [TOl shows the same critical value of \Xh \ / Pi as in Fig. [9l but for two different 
values of 60/ pi : 10~^ and 10"'^. The distance of the boundary from the line at Xjj = increases 
with 60/ Pi- 

VIII. SAWTOOTH MODELLING 

Large tokamaks can operate in regimes where the collisional fluid limit for the drift-tearing mode 
is no longer justified and the semi-coUisional theory is appropriate. The stability theory of the semi- 
collisional drift-tearing and internal kink modes is characterised by the two key parameters /3, and 
So/pi and semi-collisional effects are strong when /3 ~ 0(1) and large ion orbit effects are relevant 
when 5o/pi <^ 1. As an example, we consider sawtooth modelling in JET and ITER. For the m = \ 
mode, assuming density profiles such that n = n(0)(l — ar"^ /a?) with a ~ 1/2, we can write 

<^0 = f 77^) ^ X cm; /3 = 2/3e^— — ^— "2, where g(ri) = 1. (104) 
\ile J aq'{ri) a? [g'(ri)] 

We take the following sets of parameters 

JET : ne(0) = 10^° m^^ re(0) = 5keV, 5 = 3T, = 3m, a = Im. (105) 

ITER : ne(0) = 10^° m^^, Te{Q) = 22keV, B = 5.3T, R = 6.2m, a = 2m. (106) 
From these we derive 

JET : Pi = 0.3cm, /3e(0) = 2.2 x 10"^ v^/^ce = 3 x 10"^ (107) 

ITER : Pi = 0.4cm, ^e(O) = 3.1 x 10~^ Ve/^ce = 2 x 10"^°, (108) 

so that 

/3^^^ = 0.4/ [ag'(n)] ' , = 0.6/ [ag'(n)] ' , (109) 

5^^^ = 1.5 X IQ^^/q'in) cm, 5i^^^ = 0.4 x 10~Vg'(ri) cm. (110) 
Clearly in each case 13 ~ 0(1), while 5q/ pi <C 1, and the semi-collisional regime is appropriate^ . 

[1] When the colhsionahty is low enough to give 5o <Si de, with de = p^l \fpl the electron inertia length, the coUisionless 
regime is appropriate [2^. [2^. 
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Figure 11: Boundary of marginal stability of the dissipative internal kink mode and drift-tearing mode, 
respectively in terms of /3 and 5qA' from Ea. ([M|) (solid line) and from Eq. ([5T|) (dashed lines). The 
parameters are the same used in Fig. [HI The kink boundary does not change significantly for rje. < 2.39. A 
schematic trajectory for the evolution of a discharge following a sawtooth crash is indicated by the line with 
arrows. 



Superimposing the stability diagrams in Figs. [7] and [9] at high A'pi (where we replace —Xh in 
terms of A'r^ = —tt/Xh) as shown in Fig. [TTl we see that there are regions of parameter space 
corresponding to windows of stability to both modes. This has a potential implication for sawtooth 
modelling in the spirit of Ref. Q|, but provides more accurate stability criteria to associate with a 
trigger for the sawtooth crash. After such a crash, one expects a relatively rapid recovery, on the 
scale of the energy confinement time, of density and temperature, followed by a slower evolution of 
the q profile on the resistive diffusion time scale. An expression for A' is given by riA' = s1/6W, 
where 6W is a measure of the potential energy of the ideal internal kink mode, and si is the magnetic 
shear at ri, the radial position of the q = 1 surface. 6W is a rather complicated function of the 
pressure and current profiles and, as demonstrated in Ref. Q|, will also be a sensitive function of 
the pressure profile of fusion alpha particles in ITER. However, the contribution of the core plasma, 
as calculated in Ref. jisl for ri/a <^ 1, gives 

SW oc ( (1 - gc 
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when /3p in a measure of the poloidal f3 within the q = 1 surface. As the core plasma heats up, 5W 
reduces (due to increasing /3p) so that, provided that /3 < 0.4 (see Fig. [7]) there may be a period 
when the drift-tearing mode is unstable (a possible explanation for successor oscillations [35| after 
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a sawtooth crash), fohowed by a stable period during which the equihbrium trajectory in (/3, 5qA') 
(see Fig llip passes through a region stable to both modes. 

After the density and temperature profiles have recovered (we suppose this occurs before meeting 
the dissipative kink-instability boundary on the right of the stable region), the evolution of A' will 
be determined by that of the q{r) profile, with dependence on both si (increasing A') and ri 
(decreasing A'; however, experimental evidence suggests that the variation of ri during a sawtooth 
ramp is weak). Since 6q oc ri/si, the parameter SqA' may remain fairly constant during this phase 
if ri grows only slowly relatively to si(t). However (3 oc s^^ should decrease and the equilibrium 
trajectory could exit the stable region in (/3, 6oA') by crossing the dissipative kink boundary, as 
shown schematically in Fig[TTJ This could represent the onset of the next sawtooth collapse phase. 
The sawtooth trigger can thus be deduced from Eq. (j96p . which resembles, but is a more precise 
formulation of, the stability criterion invoked by the authors of Ref. as a prescription for a 
sawtooth trigger. 

In summary, the stability diagram in Fig. [TT] suggests a picture in which, after a sawtooth crash, 
there is first a period during which the drift-tearing/kinetic Alfven wave branches may be unstable, 
giving rise to postcursor or successor oscillations, then a phase in which these modes are stabilised 
and the dissipative internal kink is not yet unstable, until finally the threshold for the latter mode 
is reached and a sawtooth crash ensues: a repetition of this whole cycle then begins. 

The discussion above is, of course, extremely qualitative. A more meaningful technique would 
be to use a transport code to evolve the density, the temperature, the alpha particle population if 
appropriate, and the q profile during a sawtooth period and follow the resulting discharge trajectory 
relative to the stability boundaries in Fig. [TTJ Here we content ourselves with simulating the 
evolution of the q profile following a sawtooth crash, using this to calculate the evolution of /3(t). The 
simulations utilise neoclassical resistivity, corrected for finite i/,,e) the banana regime collisionality 
parameter, since this becomes important close to the magnetic axis [36|. In the case of JET, 



10-^ while for ITER, 10 ^. As an example we assume the post-crash q 

profile is given by the Kadomtsev full reconnection model js^ - other models could well have been 
invoked, such as a partial reconnection model jsl. We study a typical JET and ITER base-line 
H-mode scenario. The resulting prediction for the sawtooth periods, identified as the times that 
the instability boundary in Fig. [11] are reached as a result of the drop in /3{t), are quite plausible 
for both JET and ITER simulations: 2 sec for JET, and 200 sec for ITER. 
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IX. DISCUSSION AND CONCLUSIONS 

The stability theory of the semi-conisional drift-tearing and internal kink modes, the subject of 
this paper, is characterised by the two key parameters, /3 and 60 /pi, the latter being related to the 
collisionality parameter, C: So/pi ~ C^^'^ . Semi-collisional effects are strong when /3 ~ 0(1) and 
large ion orbit effects are relevant when 5o/pi <^ 1. Large, hot tokamaks such as JET and ITER 
operate in regimes where these conditions are satisfied. We have derived a unified theory for the 
stability of the modes at low /3 = PeiL^/ L'^)/2 and for ion Larmor orbits much larger than the 
semi-collisional layer width, i.e. pi/60 ^ 1, also extending the earlier work Q| on the kink mode. 
The effect of finite /3 on the stability of the drift-tearing and dissipative internal kink modes has 
then been elucidated and forms the basis for a sawtooth model. 

The non-local aspects resulting from the ion orbits are accommodated by working in Fourier 
transform space, where the problem reduces to solving a fourth order ordinary differential equation. 
The problem is further simplified by exploiting the two disparate scales, pi and 60, solving separately 
in the 'ion region' where the dissipative electron effects can be ignored and the 'electron region' 
where the ions can be regarded as un-magnetised. The ion region is then described by a second order 
differential equation in 'A;-space', with the ideal MHD solution involving A' providing a boundary 
condition at low k. The high-fc power solutions are then asymptotically matched to the similar 
ones emanating from the electron region solution. However, the electron region is described by a 
fourth order differential equation in A;-space and it is easier to solve this by first transforming back 
to x-space, where it becomes second order, and then calculating the Fourier transform of the small 
X asymptotic form of its solution for matching to the ion region. 

The outcome at low /3 is a unified dispersion relation that encompasses the drift-tearing mode 
and internal kink mode. The nature of the resulting modes can be usefully parameterised in terms 
of A'pif3. At low values one finds the drift-tearing mode with frequency w ~ 1 -|- cr/e, where c 
varies weakly with rje (i.e. c increases from 0.7 to 0.84 as 7]^ increases from small to large values), 
has a growth rate 7 ~ {6o/pi){A' - A'^^^^) if A' > AL^^ ~ p^'^ ^r]'^{ln{pi/6o) + c'ljl^'^), similar to 
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[see Eq. (j42l) for a precise description]. 



the stabilisation by large ion orbits shown in Ref. 
resembles that obtained for the cold ion model of Ref. 
As A'pif3 increases the drift-tearing mode is affected by a kinetic Alfven wave (KAW) propagating 
in the electron direction. Indeed for A'pif3 ~ C(l) these couple strongly and there is a cross- 
over with an exchange of stability character: the mode identified by w ~ 1 -|- crje becomes stable, 
whereas it is the KAW that becomes unstable (see Fig. [1]). Near the cross-over point the growth 
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rate is enhanced: 7 ~ {5q/ piY^"^; however, as A'pj/3 mcreases further, the growth rate diminishes 
as 7 ~ {60/ pi){A' piP^)~^] and a) — t- 1. Indeed, according to Eq. (1481) . when A' exceeds a critical 
value, A'^j-it ~ (5(^^r/e//3, it becomes stable. Stability is also ensured if i]f. > 9.4. 

We have also separately examined the stability of the w w 1 mode in the limit (3^1 and, 
somewhat surprisingly, find results closely resembling those obtained in Ref. jisl for the cold ion 
limit: the mode is completely stabilised for all A' due to shielding of the electron reconnecting 
layer as a result of the presence of radial gradients in the electron temperature. To elucidate the 
transition between instability at low /3 and stability at high j3, we have developed a treatment for 
finite /3, but one which is only strictly valid for a specific, though reasonable, value of r]e {r]e = 2.5). 
This demonstrates the onset of stabilisation due to plasma gradients when r]eP = 0.85. Below this 
critical value, /3c, the critical A' for stability is found to increase as f3 increases, as shown in Eq. ()87p . 
which is consistent with Eq. (j49p . A similar analytic calculation is also possible for the much larger 
value rje = 9.4, but in this case no unstable mode is found, irrespective of the value of f3. This 
is consistent with the threshold r/g = 9.4 predicted by low /3 theory. These stability results are 
summarised in Fig. [71 

When the ideal internal kink mode stability parameter Xh = — vr/AV^ > 0, our low /3 theory 
recovers results of Ref. Q for the effect of the large ion orbits on the ideal internal kink mode, 
although the role of the small electron dissipation is somewhat different due to the improved electron 
model used here, which allows electron temperature perturbations and hence involves the electron 
temperature gradient parameter, rje. When Xh = 0, the dissipative internal kink mode, which is 
related to the solution of our unified dispersion relation and whose stability depends entirely on 
the electron effects, remains unstable. However, we find that this mode becomes stable at a critical 
value of A', namely A'r^ = —tt/Xh = (3/ ln[{pi(3'^)/6o]. Furthermore, if rj^ > 2.38 the mode is 
stable for all A'. 

Finally we have exploited the low frequency of the dissipative internal kink mode to develop an 
analytic dispersion relation, valid for finite /3, to examine the effect of (3 on its stability. An analytic 
expression (jlOip for the marginally stable value can still be obtained. Although at first the critical 
value of —Xh reduces as f3 is increased, this trend eventually reverses and indeed all values of —Xh 
are unstable when f3 > 7/18. Figure [9] summarises these results. 

The combined stability diagram for the drift-tearing mode and the dissipative internal kink 
mode shown in Fig. [TT] provides a basis for sawtooth modelling as described in Section rVIIII In the 
aftermath of a sawtooth crash one envisages a rapid recovery of density and temperature profiles, 
followed by a slower evolution of the q profile. Essentially this leads to a decrease in the /3 parameter 



33 



until the instability boundary for the dissipative internal kink is encountered, triggering a further 
crash. This could be modelled more completely with transport codes in the spirit of Ref. ^, 
but using the more precise stability criteria given in the present theory. A simulation of the time 
evolution of the magnetic shear using neoclassical resistivity and following its consequences for the 
trajectory of /3(t) in the stability diagram, leads to plausible sawtooth periods for JET and ITER. 

At this point it is helpful to highlight our key results, expressing them in terms of the parameters 
/3 and 6q/ pi oc r^/^C^/^; the 0{1) coefficients ci — C4 below are functions of the other parameters, 
r/e, Tji, and r, although we take r = 1 in the evaluation of critical values of r/g and /3 (more detailed 
descriptions of these appear in the referenced equations): 

• At low /3 the drift-tearing mode with a; ~ 1 + 0.877e is unstable for p.jA' > ci/3r/g ln(/9j/5o) 
due to finite ion orbit effects [see Eq. (02])]. 

• At A'/9j/3 ~ 0(1) this mode interacts with the KAW and for larger values of this parameter 
it is the KAW with a; ~ 1 that is unstable, until it is stabilised when A' > C2P/I1Q [see Eq. 
(149]) ]: it is also stable for all A' if r/e > 9.4. 



For high /3 there is stability of the tearing mode for all A', the transition occurring when 
(3 = Cs/r/e [see Eq. (|86p ]: the critical A' in Eq. (|49p also increases as (3 approaches this 
transition value, as shown in Eq. (|88p . 



• The stability regimes of the drift-tearing and kinetic Alfven waves in terms of /3 and A'pj are 
summarised in Fig. [7] 

• At low 13 the dissipative internal kink mode is only unstable for A'p.;, = — 7r/>j/(rsA//) > 
C4/3/[(5o/yOjln(/32pj/(5o)] [see Eq. (J96])]; it is also stable for aU A' if ?/e > 2.4. 

• At higher /3 > 7/18, according to Eq. (jlOip there is instability for all —Xh- 

• Figures [9l and [TOl provide stability diagrams for the drift-tearing/KAW and dissipative internal 
kink modes in the space of A'(or —Xh) and f3 for higher values of A'pi. 

• Figure JTl] suggests a resulting scenario for the onset of a sawtooth crash . 

Our calculation is limited to cylindrical geometry - strictly to a sheared slab since we neglect effects 
of cylindrical curvature. In toroidal geometry one must take account of magnetic drifts. In the limit 
that the ion banana width, Wbam is small compared to 5o, one can derive similar equations to those 
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39|. While we may neglect most 



of Ref. [isl . but encompassing neoclassical transport effects 
of these for small values of the magnetic shear and the fraction of trapped particles, an important 
consequence arises from the ion neoclassical polarisation drift which reduces the parameter C by a 
factor Bp/B^ <^ 1, where Bp and Bt are the poloidal and toroidal magnetic fields. In the opposite 
limit Whan ^ ^0 one would need to include the complete drift orbits; this calculation would be 
extremely complicated but can be expected to be qualitatively similar to the large ion Larmor orbit 
case with the substitution pi — )• Whan- The fact that the effective value of C is greatly reduced in 
the torus significantly extends the experimental validity of our large orbit theory. 
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Appendix A: Internal Kink Mode 



To address the case of finite /3 for the dissipative internal kink mode, we can consider the limit 
\uj\ <ti 1, unlike the situation with the m > 1 tearing mode which has u ~ w^e- In the electron 
layer, region 2, Eq. (jlSp then becomes 



(2}/3a{s)A 



(Al) 



where 



1 + 1.71r?e + dis'^ 



(A2) 



1 + [4 - (1 + 1.71r?e)]s2 + (1 + t)ujs^ ■ 

In the region s ~ 1 we can ignore the last term in the denominator and solve Eq. ()Aip as an 
expansion in lo, retaining it only for large s ~ |ci)|~-'^/2 , when it competes with the term proportional 
to s^. Thus, setting A = Aq in lowest order we find the solution even at s =0 to be: 



A = An 



l + PCb ds'{s - s')a{s') 
Jo 



(A3) 
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At large s we find 



^ ~ Ai 



11 I— ( ci d-T , ^ 



for s —7- oo, 



witli a = 1 + l.Tlr/e, d = do — a. In the region s ~ l^^l^^^; Eq- ()Aip takes the form 

d^A _ /3 A 
dy'^ 1 + r 1 + y2 ' 

where = d^^di{l + r)a)s^, with solution 

1211 ^ 2' 4^2 2 ^/ ^^M4 2 4 2 2 ' 
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where 2-^1 is the hypergeometric function 
()A4p at small y we have 



32] and fii = a/ 1/4 + + r). Matching to solution 



61 2'^ V 



dl 1 + T 



(A7) 



To match to the ion region 1, we need the large y limit of solution (IA6 
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26ir(i-f) r(f-f) 



(A8) 



Since in the ion region we solve for the Fourier transformed current J(fc), the form of J(s) in the 
electron region can be obtained using Ampere's law. Thus, we obtain 



d 



^ 1/2 + r(-/.i) (I + S + tan [f (1 - /.i)] 

1/2 -/ii r^n) r2(|-f)s + tan[f + 



where 



1 
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U + r)(ir(f-f)r(| + f)- 



(A9) 



(AlO) 



7r(a — di/d) 

~ > 1 when \Cj\ < 1, in Eq. (jA9]) . we find that the low ^ limit of (|A9l) 



Noting that 

coincides with the low (3 calculation (??) at low a). 

Equation ()A9p . after using Eq. (j20p . must be matched to d-/d^ as derived from the region 1 
solution. An analytic expression for this ratio can be obtained by iterating the solution in |<ip/3 <^ 1, 
in a similar manner to our low /? treatment of the drift-tearing mode. Thus, in place of Eq. ([8]) we 
have 



d 
dk 



Goik) dJ 
Fo{k) dk 



(All) 
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where 



Goik) = -{l-6j-^) + Fo{k), 



(A12) 



Fo{k) = Lu-^ {ro(A;V2) - 1 - mk^^ [^k^l) - ri{e/2)] } 



(A13) 



where we retain the 0{uj) correction in Gq for large k, since in this hmit Fq ^ k ^. Solution (j22l 
is then modified to give 

du Fo{u) 



J{k) = exp 



u Go{u) 



+ Lop — / au- 



Pi Jo Gq{u) 



(A14) 



/o Go{u) Jq v^Go{v)' 
where we have expressed A' in terms of Xh- To match solution (IA14[) to that in region 2, we must first 
consider an intermediate large k range where the approximations ()A12p - ()A13p fail: k ~ ^ 1. 

In this region 



(A15) 



so that Eq. (lAlip can be written as 



^d^J dJ /3 ~ 



(A16) 



where u = ^/7^u){l + t)(1 — r]i/2) ^k, with solution 



J{k) = ai2Fi ( -/ii,-^ + - uj + a2U^ 2F1 f ^ "Mi'^ +^i;3;u ) . (A17) 



[Note, the second solution is a special case, see p. 75, Eq. (7) of Ref. 
solution ()A17p takes the form 



40|]. For \Cjk\ < 1, i.e. n < 1, 
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+ ... 
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+ a2ti^(l + ...) for \(2}k\ < 1, 

with ko = ^(l) + ^(3) — ^(3/2 — /ii)— ■;/'(3/2 + /_ii), where tp is the digamma function 
(jAlSP can be matched to the large k limit of solution ()A14p gwhere the form (jAlSP implies 
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so that 
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(in this integral we can take the — t- limit of F anf G). The result of the matching is 
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It remains to match the large u limit of solution (IA17|) 
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ai 1 + r / 

to the region 2 solution. Reintroducing k, Eq. (IA23P provides the ratio 02/01 to be matched to 
c^/c- as given by Eqs. (|A9P with the help of Eq. (|20p . i.e. 



c_ _r(-i + Mi) 
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This yields the final dispersion relation 



tan 



TT /I 



2 2+/^^ 



(A24) 



e 4' 



1 



SV^iijp, V di (1 + t)3/2 



S(a;)+tan[f (^-/ii)] 
H(a;)+tan[f (i + /ii)] 



1/2 -//I r=^(pi) S(^) + z^e 
l/2 + /.ir2(-^0s(^) + ,4-e- 



where T,{uj) = 202/01. We note that we can simplify Eq. ()A25P since |H(a))| oc 
by Eq. ^MOil, 



-1/2 



(A25) 



^ 1, as given 



80FwpiV di (1 + t)3/2 



2mi 



1/2 -pi r2(pi) S(^)+i4^e 



(A26) 



1/2 + pi r2(-pi)2(^) + i^e- 

Equation ()A26P provides a finite /3 dispersion relation for the dissipative internal kink mode, where 
we recall 5o/pi <^ 1 and pi is independent of u. 
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